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ABSTRACT 

To empirically calibrate the IR surface brightness fluctuation (SBF) distance 
scale and probe the properties of unresolved stellar populations, we measured 
fluctuations in 65 galaxies using NICMOS on the Hubble Space Telescope. The 
early-type galaxies in this sample include elliptical and SO galaxies and spiral 
bulges in a variety of environments. Absolute fluctuation magnitudes in the 
F160W (1.6 /xm) filter (Mpieow) were derived for each galaxy using previously- 
measured /-band SBF and Cepheid variable star distances. F160W SBFs can be 
used to measure distances to early- type galaxies with a relative accuracy of ~10% 
provided that the galaxy color is known to ~0.035 mag or better. Near-IR fluc- 
tuations can also reveal the properties of the most luminous stellar populations 
in galaxies. Comparison of F160W fluctuation magnitudes and optical colors to 
stellar population model predictions suggests that bluer elliptical and SO galaxies 
have significantly younger populations than redder ones, and may also be more 
metal-rich. There arc no galaxies in this sample with fiuctuation magnitudes con- 
sistent with old, metal-poor (t > 5 Gyr, [Fe/H] < —0.7) stellar population models. 
Composite stellar population models imply that bright fluctuations in the bluer 
galaxies may be the result of an episode of recent star formation in a fraction of 
the total mass of a galaxy. Age estimates from the F160W fluctuation magni- 
tudes are consistent with those measured using the H/? Balmer line index. The 
two types of measurements make use of completely different techniques and are 
sensitive to stars in different evolutionary phases. Both techniques reveal the 
presence of intermediate-age stars in the early-type galaxies of this sample. 

Subject headings: galaxies: evolution — galaxies: stellar content — galaxies: 
distances and redshifts 



1. Introduction 

The techniques for measuring surface brightness ffuctuations (SBFs) were developed 
primarily with the goal of determining extragalactic distances. The optical SBF method 
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has proven to be remarkably useful for estimating distances (see Blakeslee et al. 1999 for a 
review). The initial motivation for extending ground-based SBF measurements to near-IR 
wavelengths was to take advantage of intrinsically brighter fluctuations and better atmo- 
spheric seeing to reach much greater distances (Jensen et al. 2001, hereafter J2001; Jensen, 
Tonry, & Luppino 1999; Liu & Graham 2001). With Hubble Space Telescope's (HST) Near 
Infrared Camera and Multi-Object Spectrometer (NICMOS), the low background and lack 
of atmospheric seeing offer significant advantages over ground-based IR SBF measurements. 

Accurate SBF distance measurements rely on the empirical calibration of absolute fluc- 
tuation amplitudes. Several studies have used observations of galaxies with previously- 
measured distances to determine absolute fluctuation magnitudes, which were found to agree 
with the predictions of stellar population models. The first IR SBF measurements were made 
by Luppino & Tonry (1993), who found that M32 has a somewhat brighter (~0.25 mag) K- 
band SBF magnitude than the bulge of M31. Pahre & Mould (1994) measured X-band SBF 
magnitudes for a sample of Virgo galaxies, and found that two of eight had significantly 
brighter (>1 mag) fluctuations than the others. Pahre & Mould's results were confirmed by 
Jensen, Luppino, & Tonry (1996) in their sample of seven Virgo cluster elliptical galaxies. 
The IR SBF signal for one of the two anomalous galaxies (NGC 4365) was later found to 
have been over-estimated due to the contribution of undetected globular clusters (Jensen, 
Tonry, & Luppino 1998); the other (NGC 4489) still appeared to have brighter fiuctuations 
than the others. 

To the 1996 Virgo sample, Jensen et al. (1998) added galaxies in the Fornax and 
Eridanus clusters. Three of the bluest galaxies in their sample showed i^-band fluctuation 
magnitudes that were ~0.25 mag brighter than the others. Jensen et al. noted that the 
models implied younger stellar populations in these bluer ellipticals. However, the sample 
size was small and the range in (V—I) color was limited; hence, the slope they measured 
was not statistically significant and they adopted a constant K-hand SBF calibration. The 
K-hand fluctuations for the bluest galaxy in their sample (NGC 4489) were signiflcantly 
brighter than the others in the Virgo cluster, but the measurement was not trusted due to 
its low signal-to-noise (S/N) ratio. Mei, Silva, & Quinn (2001c) subsequently re-observed 
this galaxy and found comparably bright A'-band fluctuations. 

Recently, a larger sample of K-hand SBF magnitudes for Fornax cluster galaxies has 
been measured by Liu, Graham, & Chariot (2002). Combined with the previously published 
data, the resulting sample covered a much larger range of galaxy properties (e.g., color and 
luminosity). Liu et al. (2002) discovered that the brighter i^T-band fluctuations seen in bluer 
galaxies in their sample and in earlier data sets (Luppino & Tonry 1993; Pahre Sz Mould 
1994; Jensen et al. 1998) were correlated with (V—I). Liu et al. found the X-band slope 
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with color to be comparable to that for /-band SBFs. This relation naturally explained 
the bright IR SBFs previously measured in M32 (Luppino & Tonry 1993) and NGC 4489 
(Mei et al. 2001c). Moreover, a generic prediction of IR SEP models is that IR fluctuations 
are sensitive to variations in the ages and metallicities of stellar populations. Hence, the 
discovery of a systematic relation between IR SBFs and galaxy color means that age and 
metallicity are related in a way that reflects the star formation history of a galaxy. Liu et 
al. (2002) concluded that early-type galaxies in clusters have a significant spread in age and 
approximately solar metallicities. 

Liu et al. (2002) also reported finding three Fornax cluster galaxies with i^'-band fiuctu- 
ation magnitudes brighter than the other galaxies with similar (V—I) colors. These galaxies 
were not included in their calibration fit. They are primarily lower luminosity galaxies, and 
models suggest the presence of a high-metallicity burst of star formation in the last few Gyr. 
Mei et al. (2001a) have confirmed that one of the three Fornax galaxies (NGC 1427) has 
unusually bright K-hand fiuctuations. 

To calibrate the F160W (1.6 /xm) SBF distance scale and better understand the nature 
of the bluer, low-luminosity eUiptical galaxies, we measured F160W SBF magnitudes in a 
large sample of 65 galaxies spanning a wide range in color {{V—I)o — 1.05 to 1.28). All 
the data presented here were obtained using the NIC2 camera of NICMOS. The data are 
of uniform image quality, and the S/N ratios are large. Many of the galaxy images were 
collected for other programs and retrieved from the public archive. 

Distances to the galaxies in this sample were taken from the growing collection of 

Cepheid variable star distances measured with the HST (Freedman et al. 2001, Ferrarese 
et al. 2000b; Gibson & Stetson 2001; Saha et al. 2001), and from the extensive /-band SBF 
survey (Tonry et al. 2001; Ajhar et al. 1997; Lauer et al. 1998; Tonry et al. 1997). The /- 
band SBF distances used here were calibrated using Cepheid distances to the same galaxies; 
no assumptions about group or cluster membership were made (Tonry et al. 2001). For the 
current study we shifted the /-SBF distances published by Tonry et al. by —0.16 mag to the 
new Cepheid zero point of Freedman et al. (2001), which makes use of the improved Cepheid 
period- luminosity relations published by Udalski et al. (1999). The distances presented in 
this paper, whether from Cepheids directly or from /-band SBFs, are subject to the sys- 
tematic uncertainties in the Cepheid distance scale (Ferrarese et al. 2000a). We adopted a 
distance modulus to the Large Magellanic Cloud (LMC) of 18.50 mag (Freedman et al. 2001; 
Carretta et al. 2000), and acknowledge that the continuing debate over the LMC distance 
remains one of the largest sources of systematic uncertainty. 

In this paper we present a new calibration of the F160W distance scale and explore its 
sensitivity to galaxy color. We also compare the absolute fiuctuation magnitudes to those 
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predicted by three sets of stellar population models: the widely-used models of Worthey 
(1994), the Bruzual & Chariot (1993, 2002) models (as published by Liu, Chariot, & Graham 
2000 and Liu et al. 2002), and the Vazdekis (1999, 2001) models (as published by Blakeslee, 
Vazdekis, & Ajhar 2001). SBF ages determined using stellar population models are compared 
with those measured using the age-sensitive H/3 Balmer line strength. 

2. Observations 

2.1. F160W SBF Measurements 

The observations of galaxies used in this study come from two basic types of NICMOS 
program, each of which accounts for roughly half the data. The first set of observations came 
from two programs explicitly designed for SBF measurements. The first of these (NICMOS 
program ID 7453) targeted galaxies in Virgo, Leo, and Fornax to empirically cahbrate the 
F160W SBF distance scale for more distant measurements. The second (program ID 7458) 
measured fluctuations in a large sample of Fornax galaxies at multiple wavelengths to better 
understand stellar populations. A subset of the measurements presented in this paper was 
published by J2001. The second category of observations is comprised of data taken from 
the public archive from a variety of programs. Most of these are short, "snapshot" survey 
images of the centers of galaxies, and have only a minimal number of individual exposures. 
Observational data are hsted in Table 1. 

The observations presented in this study were taken with the N1C2 camera (19.2 arcsec 
fleld of view) through the F160W filter. F160W N1C2 observations of approximately 300 
galaxies with heliocentric redshifts less than 10,000 kms^^ were retrieved from the public 
archive and examined. A smooth fit to each galaxy was subtracted from the pipeline-reduced 
image and the residual frame examined. All the galaxies that were judged to be hopelessly 
dusty were rejected from further consideration. Data for the remaining ~200 galaxies were 
reduced again from the raw data and the SBF analysis completed following the same proce- 
dures as described by J2001. The subset of 65 galaxies presented here arc those for which 
reliable distances are known from either Cepheid variable stars or /-band SBFs. 

The methodology for determining F160W fiuctuation magnitudes was very similar to 
that described by J2001, with some simplifications. We used the software developed by the 
NICMOS GTO team to prepare the images for analysis (Thompson et al. 1999). Dark 
current was first subtracted from the raw N1C2 F160W images. The multiple reads of each 
MULTIACCUM sequence were combined and cosmic rays identified and removed. Adjust- 
ment of the bias pedestal for each quadrant was performed as described by J2001, although 
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it was unnecessary in most cases. Exposures from individual telescope pointings were flat 
field corrected and combined. Residual images from cosmic rays were not a significant source 
of contamination in these high-S/N ratio data. 

The photometric zero point used by J2001 was determined by M. Rieke and the NICMOS 
team. Additional standard star measurements combined with a better tie to ground-based 
photometry has yielded a better zero point for the NIC2 F160W filter. The new photometric 
zero point that we adopted for this study is 0.033 mag fainter than that used by J2001. The 
new calibration for N1C2 is 2.126x10^^ JyADU^^s^^. The magnitude zero point on the 
Vega system is 1083 Jy. This new calibration is within the uncertainty in zero point published 
by J2001. If the current zero point were applied to the J2001 data without any other changes 
to the calibration, the resulting fluctuation magnitudes would be 0.03 mag fainter, and the 
Hubble constant would be larger by 1.6%. SBF magnitudes were corrected for Galactic 
extinction using the measurements of Schlegel, Finkbeiner, & Davis (1998). We adopted 
Ab^4:M5E{B-V) and Ah = 0.132^b (Schlegel et al 1998). 

The galaxies in this sample are much closer than the distant galaxies used by J2001 to 
determine the Hubble constant (5 included here were taken from the intermediate-distance 
sample of J2001). The median S/N ratio was 16 per pixel for this data set, which, given 
the large number of pixels sampled, is more than sufficient to achieve a highly-reliable mea- 
surement. Most of the data had S/N ratios between 10 and 20; the full range includes 
measurements with S/N ratios as low as 5 and as high as 4000. The fluctuation power 
was determined by fitting the scaled power spectrum of the reference point-spread function 
(PSF) to the spatial power spectrum of the cleaned and galaxy-subtracted data. The library 
of empirical PSF stars collected as part of the IR SBF Hubble constant project (J2001) 
was used to perform the SBF analysis. A full discussion of the techniques for calculating 
fluctuation magnitudes mpieow from NICMOS data is presented by J2001. 

We found that correcting the F160W SBF magnitudes for undetected globular clusters or 
background galaxies was unnecessary for the relatively-nearby galaxies of this study because 
the stellar SBF signal always dominated over other sources of variance. The flnal galaxy- 
subtracted images were masked of visible point sources and dusty regions before proceeding 
with the SBF analysis. Only the five intermediate-distance galaxies from J2001 required 
corrections for undetected globular clusters, background galaxies, and residual background 
patterns (as described by J2001). 

The uncertainties in apparent fluctuation magnitudes mpieow were typically 0.1 mag 
or less. The primary components of the uncertainty were the fit to fiuctuation power and 
the PSF normalization (usually 0.05 mag each). The contribution from sky subtraction 
was much smaller (typically 0.01 mag). Contributions to the uncertainty from undetected 
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globular clusters or background galaxies were negligible (<0.01 mag), and contaminating 
power from residual cosmic rays or incomplete bias subtraction were unmeasurable in these 
high-S/N ratio measurements. 

2.2. Distances and Absolute Fluctuation Magnitudes 

Computing the absolute fluctuation magnitudes Mpieow for the galaxies in our sample 
required independent distance measurements. All the distance moduli used in this study 
(Table 2) were based, either directly or indirectly, on Cepheid variable star distances. Two 
classes of distance measurements are presented here: first, a large set of 61 /-band SBF 
distances (calibrated using Cepheid distances as described by Tonry et al. 1997), and second, 
a smaller set of nine Cepheid distances measured using WFPC2 on the HST. 

Because only a few of the spiral galaxies with measured Cepheid distances have smooth, 
dust-free regions appropriate for SBF analysis, the majority of the absolute fluctuation mag- 
nitudes presented in this paper were computed using /-band SBF distances (Tonry et al. 
2001). Most of the /-band SBF distance measurements were made using ground-based tele- 
scopes, but five galaxies were observed using WFPC2 (J2001; Lauer et al. 1998; Ajhar et 
al. 1997). The empirical /-band SBF cahbration (Tonry et al. 1997) adopted for this study 
used Cepheid and /-SBF distance measurements to seven galaxies for which both types of 
measurement are possible. No assumptions about group or cluster membership were made to 
connect the Cepheid and /-band SBF distance scales. The original Key Project calibration 
of the Cepheid distance scale (Ferrarese et al. 2000b, Freedman & Madore 1990) was used for 
the empirical /-band SBF calibration presented by Tonry et al. (1997, 2001). For the /-band 
SBF distances presented here we have updated the zero point using the greatly-improved 
period-luminosity relations determined by Udalski et al. (1999) using 650 LMC Cepheids. 
The new period-luminosity relation results in a shift of —0.16 mag in the distance moduli of 
all the /-band SBF galaxies. No metaUicity correction to the Cepheid distances is adopted 
for this study. 

The modifications to the period-luminosity relation adopted by Freedman et al. (2001) 
are distance dependent. The new period-luminosity relation is different for the /-band, while 
at V it is unchanged (Udalski et al. 1999). The effect is that the new period-luminosity 
relation predicts higher reddening corrections for redder Cepheids, and therefore smaller 
distances. The effect is largest in Cepheids with longer periods. Because only the brightest 
(longest-period) Cepheids are detected in the most distant galaxies, the offset between the 
previously-published distances and those derived using the new calibration are largest in 
the most distant galaxies (~0.2 mag, or 10% in distance). The offset to the /-band SBF 
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calibration due to the new period- luminosity relation is 0.16 mag, or 8% in distance. 

Nine galaxies in the F160W NICMOS sample have Cepheid distances measured using 
HST (Ferrarese et al. 2000b; Preedman et al. 1994, 2001; Gibson & Stetson 2001; Gibson et 
al. 1999, 2000; Graham et al. 1997; Saha et al. 1996, 2001). A distance modulus to the LMC 
of 18.50 mag was adopted for the current study and for the Key Project papers (Ferrarese et 
al. 2000b; Freedman et al. 2001). This distance to the LMC appears justified in light of many 
recent distance measurements (e.g., Carretta et al. 2000), although significant differences 
certainly remain between techniques and investigators (Udalski et al. 1999; others). A full 
discussion of the LMC distance issue is beyond the scope of this paper; future adjustments to 
the LMC distance modulus of 18.50 should be applied as a constant offset to all the distances 
used herein and to the resulting F160W SBF calibration. The conclusions of this paper are 
not affected by an uncertainty of 0.15 mag in the LMC distance. 

Freedman et al. (2001) also applied an empirical metallicity correction of —0.2 ±0.2 
mag dex~^ in metallicity to the Cepheid calibration such that more metal-rich Cepheids 
are intrinsically brighter (Freedman et al. 2001; Gibson & Stetson 2001; Kennicutt et al. 
1998). For comparison, separate distances for the Key Project galaxies are listed in Table 2 
using the new period-luminosity relation alone ("new PL", Freedman et al. 2001; Gibson & 
Stetson 2001), and using both the new period-luminosity relation and empirical metallicity 
correction ("new PL+Z", Freedman et al. 2001). We chose not to adopt the empirical metal- 
licity correction to the Cepheid distances endorsed by Freedman et al. (2001) because it is 
not yet entirely clear that the metallicity correction is justified. Udalski et al. (2001) find 
no evidence for a trend in the luminosity of Cepheids with metallicity. Another study based 
on theoretical models of Cepheid structure predicts a correction to the period-luminosity 
relation of -1-0.27 mag dex~^ in metallicity, which is similar in magnitude but opposite in 
sign from the empirical relation adopted by Freedman et al. (Caputo, Marconi, & Musella 
2002). The apparently better agreement with the reliable maser distance to the galaxy NGC 
4258 is offered as evidence that the theoretical relationship is more realistic. Of course a 
slightly smaller distance modulus to the LMC would also explain the difference, so the maser 
distance to NGC 4258 cannot be regarded as evidence that the metallicity correction must 
have a particular sign or magnitude. Given the uncertainty in the metallicity correction at 
this point, we chose not to adjust the /-band SBF calibration. In the end, the size of the 
metallicity correction is not very important provided that the correction is applied consis- 
tently through the LMC rung of the distance ladder. We note that the metallicity-corrected 
Cepheid distances reported by Freedman et al. (2001) are based on Cepheid measurements 
of LMC Cepheids that have not been corrected for metallicity. Given the similarity in metal- 
licity between Galactic Cepheids and most of the spiral galaxies in which Cepheids have been 
detected, the effect of a fully-consistent metallicity-corrected Cepheid distance scale and the 
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uncorrected calibration we adopted for this paper is minimal (approximately 0.02 mag). The 
data plotted in the figures are based on distances that include no metallicity corrections. 

The uncertainties in mpieow were combined in quadrature with the uncertainties in dis- 
tance modulus to get the final uncertainties in Mpieow presented in Table 2. The median 
uncertainty in Mpieow was 0.18 mag for this sample. The estimated systematic error in the 
Cepheid distance scale, which wc incur regardless of which distances are used to compute 
MpieoW) is 0.16 mag (Ferrarese et al. 2000a), but could be larger due to blending or other 
effects (Gibson & Stetson 2001; Mochejska et al. 2000; Ferrarese et al. 2000c). The dom- 
inant sources of systematic uncertainty are the distance to the LMC (0.13 mag) and the 
WFPC2 photometric zero point (0.09 mag). The systematic uncertainty is not included in 
the uncertainty in Mpieow listed in Table 2. It is important to note that the uncertainty 
in Mpieow is correlated with the uncertainty in (V—I) for the values derived from 7-band 
SBF distances. Mpieow is a function of the 7-SBF distance modulus, which is a function of 
(V-I) color (Tonry et al. 1997). 

2.3. Galaxy (V-I) Colors 

To calibrate the F160W SBF distance scale, we used the optical {V—I)o color corrected 
for extinction to constrain the fiuctuation magnitude dependence on stellar population. Most 
of the {V—I)o color data in Table 2 were taken from the optical /-band SBF survey (Tonry 
et al. 2001). Galactic extinction corrections were made using the Schlegel et al. (1998) 
extinction maps. The data were collected using ground-based telescopes in annular regions 
that were typically much larger than the central 20-arcsec regions imaged using N1C2. In 
a few cases we checked the optical colors within the N1C2 field of view to ensure that color 
gradients within individual galaxies were not significant. For the early- type galaxies in this 
sample, color gradients do not produce significant color differences between the F160W and 
/-band regions. A few galaxies, including the spirals not included in the /-band SBF survey 
(Tonry et al. 2001), were re-observed to allow a direct {V—I)o color measurement within 
the NIG2 field of view. 



3. Using F160W SBF Magnitudes to Measure Distances 

The primary motivation for exploring the variation of Mpieow with stellar population 
was to better cahbrate F160W SBFs as a distance indicator. With the high spatial resolution 
and low background achieved using NIGMOS, J2001 demonstrated that fiuctuations could 
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be measured in modest exposures (one or two orbits) to distances beyond 100 Mpc. J2001 
utilized a limited subset of the current data set and the Cepheid zero point calibration 
of Ferrarese et al. (2000b) to empirically calibrate the F160W SBF distance scale. They 
measured distances to a sample of galaxies reaching redshifts of 10,000 kms"^ for the purpose 
of determining the Hubble constant. J2001 restricted their calibration to galaxies redder than 
{V—I)o > 1.16 and found no significant slope of Mpieow with galaxy color. 

A constant Mpieow calibration is inappropriate for the full color range spanned by the 
galaxies in this sample (Fig. 1). For the full data set, wc fitted a slope adopting the maximum 
likelihood method described by Liu ct al. (2002) to account for the correlated uncertainties 
in Mpieow and {V—I). Uncertainties are correlated because Mpieow depends on /-band SBF 
distances, which are in turn computed using the galaxy {V—I) color. The Liu et al. (2002) 
method accounts for this non-zero covariance between Mpieow a-^d (V—I) in determining 
the best-fitting slope and intercept. If the covariance were ignored, the slope of the fitted 
fine would be biased (see Liu et al. 2002 for a discussion). 

We chose to restrict the sample of galaxies used for cafibration purposes to those that 
show no sign of dust in the NIC2 field of view because clumpy dust makes a galaxy look 
bumpier, and hence fiuctuation magnitudes brighter. Many of the rejected galaxies have 
well-defined dust lanes that can be masked. The SBF magnitudes are most likely unaffected 
by the dust, but we exclude them from the calibration fit to be safe. Given the large size of 
the sample, we can afford to be conservative in cxchiding dusty galaxies for the purpose of 
computing the calibration. Fitting the 47 dust-free galaxies gives 

^Fi60W = (-4.86 ±0.03) + (5.1 ±0.5) [(\/-/)o- 1.16] for 1.05 < (\/-/)o < 1.24. (1) 

The per degree of freedom is 1.17, indicating that the uncertainties are reasonable. The 
slope of 5.1 ±0.5 mag per mag in (V—I) color is steep, and similar to the slope of 4.5 ±0.25 
measured at /-band (Tonry et al. 1997). Liu et al. (2002) found a slope of 3.6 ±0.8 at 
K (2.2 fim). Stellar population models initially predicted that the slope of IR fiuctuation 
magnitude with color would be opposite that at /-band (Worthey 1994), but it is now clear 
that a significant slope of the same sign persists to 2.5 /im. The implications of a positive 
slope are discussed in Section 4. 

The relation in Equation 1 is based on /-band SBF distances calibrated using the im- 
proved Cepheid period-luminosity relation of Udalski ct al. (1999). The Cepheid distances 
from Freedman et al. (2001) without metallicity corrections were adopted, and we used 
the same LMC distance modulus of 18.50 mag as Freedman et al. In Figure 1 symbols of 
three sizes are plotted, with the largest symbols for galaxies with uncertainties less than 0.2 
mag and the smallest for measurements with uncertainties greater than 0.3 mag. Median 
uncertainties in each bin are plotted. Three galaxies in the "dust-free" sub-sample also have 
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Cepheid variable star distances. They were not included in the fit because they are not 
independent of the /-band SBF measurements, and their uncertainties arise from different 
sources. They are, however, plotted in Figure 1 to show that the calibration derived using 
/-band SBF distances agrees with direct Cepheid distance measurements. 

The value of {V—I)o — 1.16 marks the color at which the strong slope observed in the 
bluer galaxies makes way to an apparent flattening at the red end. Fits restricted to galaxies 
with {V—I)o> 1.16 are consistent with a slope of zero, as found by J2001 using a smaller 
calibration sample. The sloping fit over the entire color range spanned by the sample, 
however, is statistically robust. The brightest cluster galaxies observed by J2001 are all 
redder than {V—I)o = 1.16, and the use of a constant Mpieow by J2001 was justified. The 
stellar population models presented in the following sections suggest that the steep slope 
shown in Figure 1 cannot continue to arbitrarily red values of (y— /)because galaxies cannot 
be made of arbitrarily old stars. We therefore adopt the relation between M-pimw and 
{V—I)q given in equation 1 and emphasize that it is only applicable to galaxies with (V—I) 
colors between 1.05 and 1.24. 

The small statistical uncertainty in the zero point and the value of close to 1 suggests 
that the "cosmic scatter", or the variation in Mpieow not accounted for with the single 
(V—I) parameter, is small. The calibration presented here suggests that relative distances to 
galaxies with (V—I) colors typical of elliptical galaxies can be measured with 10% accuracy 
or better provided that the uncertainty in {V~I)o is 0.035 mag or smaller, and that the 
photometric zero point is known to 0.05 mag or better. This estimate of the statistical 
uncertainty does not include the systematic uncertainty in the distance scale zero point, 
which includes the uncertainties in the Cepheid distance scale and the distance to the LMC. 
The scatter increases significantly when dusty galaxies are included. The scatter below the 
best-fit line in Figure 1 appears larger than above it. This could be caused by the F160W SBF 
amplitude being under-estimated, or the /-band SBF amplitude being over-estimated (and 
the distances under-estimated). The uncertainties in the /-band SBF magnitudes are larger 
than 0.2 mag for most of the outlying points, but there is no obvious sign of a systematic 
problem with the /-band SBF measurements. The fit has been appropriately weighted by 
the uncertainties, and is consistent with the most accurate measurements. Excluding them 
would not significantly change the relation in Equation 1. 

J2001 used a constant fiuctuation magnitude Mpieow = —4.86 in their determination of 
the Hubble constant. It is fair to ask what effect the new sloping calibration and revised 
Cepheid calibration would have if they had been adopted instead. To address this issue, 
we selected the five galaxies from the J2001 sample that fall within the valid color range 
{V—I)o < 1.24 of the current distance calibration and calculated the Hubble constant for the 
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preferred flow models using exactly the same techniques used by J2001. The increase in the 
Hubble constant that results from the changes in the zero point and slope of the calibration 
is not sensitive to the exact flow model adopted or rejection of low-5'/iV ratio observations 
(sec J2001 for a discussion of these effects). Because the new Cepheid calibration and NIC2 
zero point result in a significantly fainter calibration zero point, the value of the Hubble 
constant would increase by 10% from 76 or 77 to 85 kms^^Mpc"^. If the new calibration 
is applied to only the most distant galaxies from J2001, the Hubble constant would increase 
from 72 kms~^ Mpc~^ to 79. Only a small fraction (1.5% of the 10%) of the change is due 
to the slope of the calibration; the 0.16 mag change in the Cepheid zero point resulting 
from the adoption of the improved period-luminosity relation is the dominant factor. If we 
had adopted the mctallicity corrections endorsed by Preedman et al. (2001), the increase in 
the Hubble constant would have been only 5% (corresponding to a zero point shift of 0.06 
mag instead of 0.16 mag). All the changes described happen to affect the calibration in the 
same direction, i.e., reducing the absolute brightness of the fluctuations and decreasing the 
implied distances. 



4. Stellar Population Models 

Infrared SBF magnitudes provide important new constraints on theoretical and semi- 
empirical stellar population models. The models provide insight into why fluctuation mag- 
nitudes vary with color the way they do. Each single-burst, constant-metallicity model is 
constructed by adopting an initial mass function and uniform composition for an ensemble 
of theoretical stars. These stars are then allowed to evolve according to the constraints 
and assumptions of the particular model. At each time step, the luminosity function of the 
ensemble of stars is integrated to determine observable values, including broad-band colors, 
line indices, and fluctuation magnitudes. Luminosity fluctuations in a galaxy are propor- 
tional to the second moment of the luminosity function (Tonry & Schneider, 1988), and can 
be computed from the model luminosity function at each time step. Because fluctuations 
are weighted to the most luminous stars in the population, SBFs provide a way to better 
measure the contributions of luminous first-ascent red giants and asymptotic giant branch 
(AGE) stars in unresolved stellar populations. 

The first detailed models to successfully predict SBF magnitudes were those of Worthey 
(1994). SBF predictions from the models of Bruzual & Chariot (Liu ct al. 2000, 2002) and 
Vazdckis (Blakeslee et al. 2001) have recently been published. In this paper, we compare 
empirical F160W fluctuation magnitudes Mpieow to the three sets of models to learn more 
about the relative ages and metallicities of the galaxies in our sample. These three sets of 
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Fig. 1. — Absolute fluctuation magnitudes Mpieow plotted vs. the extinction-corrected 
(y— 7)o color for the 47 galaxies that show no signs of dust in the NIC2 field of view. The 
square points represent Mpieow values derived using ground-based /-band SBF distances, 
and the triangles indicate galaxies with J-band SBF distances measured with WFPC2 on 
the HST. The size of each point indicates the uncertainty in Mpieow- The largest points 
have uncertainties less than 0.2 mag, the medium-sized points fall between 0.2 and 0.3 mag, 
and the smallest points have uncertainties greater than 0.3 mag. Median error bars for each 
point size are shown at the top of the figure. The circles indicate three galaxies with reliable 
Cepheid distances (NGC 224=M31, NGC 3031=M81, and NGC 4725). They are also plotted 
using their 7-band SBF distances. 
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models are the most successful models in common use today for which SBF magnitudes have 
been published. We also compare the model predictions to the measured fluctuation color 

(m/— mFi60w)j which is independent of the distance measurements. Agreement between 
observed fluctuation magnitudes and stellar population model predictions can therefore be 
used to simultaneously conflrm the reliability of SBF distance measurements and reinforce 
the age and metallicity interpretations. 

4.1. The Bruzual &: Chariot Models 

Absolute F160W fluctuation magnitudes are plotted as a function of extinction-corrected 
{V—I)o color in the top panel of Figure 2. Data for all 65 galaxies are plotted, including 
galaxies containing dust that were excluded from the calibration flts. Points for galaxies 
with Cepheid distances are included as well. In the lower panel of Figure 2 we plotted 
the fluctuation color (mi — mpieow) as a function of galaxy {V—I) color. The fluctuation 
color is independent of distance, and therefore provides additional insight into the stellar 
populations of the galaxies without being subject to the systematic uncertainties of the 
distance measurements. 

Theoretical fluctuation magnitudes and colors from Liu et al. (2000, 2002) are shown 
in Figure 2 as dashed (constant age) and dotted (constant metallicity) hues. The Liu et 
al. (2002) SBF predictions were based on the single-burst population synthesis models of 
Bruzual & Chariot (1993, 2002), originally described by Liu et al. (2000), and then slightly 
improved by Liu et al. (2002). The models provide several choices for the input evolutionary 
tracks and spectral libraries; a detailed comparison of the various options is given by Liu et 
al. (2000). We use their preferred set of inputs: the Padova evolutionary tracks of Bertelli 
et al. (1994), semi-empirical spectral energy distributions from Lejeune, Cuisinier, & Buser 
(1997), and a Salpeter initial mass function. 

The Bruzual & Chariot models (Liu et al. 2002) imply that the trend towards brighter 
fluctuations in bluer galaxies is consistent with younger stellar population models with high 
metallicities. The redder galaxies appear older and more metal-poor than the bluer ones 
(model ages are relative to the age of the Universe, thus a model age of 17 Gyr does not 
imply that a population formed before the Universe was born). A similar trend has been 
seen in X-band fluctuations (Jensen et al. 1998; Mei, Quinn, & Silva 2001b; Liu et al. 
2002). It is interesting to note that there are no fluctuation magnitudes among the bluer 
galaxies that are consistent with stellar population models older than about 5 Gyr. None 
of the galaxies appear to have metallicities lower than approximately [Fe/H] = —0.7. The 
spread in fluctuation magnitudes with color traces a line of models with nearly constant 
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metallicity down to (y— /)o = 1.15; bluer galaxies appear more metal-rich than those with 
(y-7)o>1.15. 

The distance-independent plot of fluctuation colors in the lower panel of Figure 2 leads 
to the same conclusion about the relative ages of galaxies with different colors. The absolute 
metallicities implied by the (m/ — mpieow) niodels are very similar to those in the upper panel. 
The bluer galaxies have the same fluctuation magnitudes as stellar population models with 
significantly younger ages and slightly higher metallicities than redder galaxies. Agreement 
between the Mpieow and the distance-independent (m/ — mpieow) niodcls suggests that the 
total systematic error in the distance scale calibration is likely to be of order 0.1 mag. 



4.2. The Vazdekis Models 

In Figure 3 wc compare our data to the F160W SBF magnitudes and {V—I)o colors 
computed for the Vazdekis models (Blakcslcc ct al. 2001; Vazdekis 2001). Additional mod- 
els for populations younger than 4 Gyr were retrieved from http://www.iac.es/galeria 
/vazdekis/col_lick.html. The Blakeslcc ct al. models make use of the new Padova 
isochrones of Girardi et al. (2000), which are transformed to the observable plane using the 
empirical stellar libraries of Alonso, Arribas, & Martinez-Roger (1996, 1999) and Lejeune et 
al. (1997, 1998). The Blakeslee et al. (2001) SBF models were computed for the //-filter, 
so it is necessary to shift the models to F160W. We adopted an empirical correction to the 
absolute fluctuation magnitudes of 

Mfi60W = Mh + 0.10(Mj-Mk) (2) 

based on photometry of red main sequence stars (Stephens et al. 2000). This correction, 
which is typically 0.17 mag for most models, has the effect of shifting the model lines in 
Figure 3 down with respect to the points. We also computed the offset between H and 
F160W by convolving the filter profiles with synthetic spectra of red giant stars (B. Plez, 
private communication; see also Bessell, CasteUi, & Plez 1998). The resulting relation 

Mpieow = Mh + 0.116(M j-Mk) + 0.026 (3) 

is very close to that derived from the Stephens et al. (2000) data and has no effect on the 
interpretation, so we adopted the empirical correction in Equation 2. 

The top panel in Figure 3 reveals very good agreement between the models as corrected 
using Equation 2 and the data calibrated using the distances based on the Freedman et 
al. (2001) Cepheid cahbration without metaUicity corrections. Uncertainties in the distance 
calibration are avoided when we compare the fluctuation color (m/ — mpieow) fo the models. 
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Fig. 2. — The upper panel shows absolute fluctuation magnitudes plotted vs. galaxy optical 
color. The circular points were calculated using Cepheid distances, the square points using 
ground-based /-band SBF distances, and the triangles using HST 7-band SBF distances. 
Open symbols indicate galaxies that show signs of dust. Clumpy dust can cause fluctuations 
to appear brighter. The lines behind the points indicate Liu ct al. (2000, 2002) models of 
constant age (dashed lines) and metallicity (dotted lines). The ages (Gyr) and mctallici- 
ties ([Fc/H]) of the models arc indicated. The lower panel shows the distance-independent 
fluctuation color as a function of (V^— /)o. 
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The ages and metallicities implied by the {mj — mpieow) models are in very good agreement 
with those of the Mpieow models and with the Liu et al. (2002) models. Once again, a change 
in the distance scale zero point of more than approximately 0.1 mag would compromise the 
compatibility between the models. 

The Vazdekis models show a compression in Mpieow of the higher-metaUicity and 
younger-age tracks compared to the Liu et al. (2000, 2002) models. This is partly due 
to the fact that the maximum metallicity shown for the Blakeslee et al. (2001) models is 
0.2, while that of Liu et al. is 0.4. It is also partly due to the difference in slope between 
the youngest populations. While there are many differences in the details of the Vazdekis 
and Bruzual & Chariot models, the choice of isochrones is perhaps the most significant for 
the current study. If we apply the newer evolutionary tracks (Girardi et al. 2000) adopted 
by Blakeslee et al. (2001) to the Bruzual & Chariot models (Liu et al. 2000, 2002), we find 
that the resulting age tracks are nearly horizontal and lie midway between the two sets of 
models shown in Figures 2 and 3. 



4.3. The Worthey Models 

Figure 4 shows the F160W data plotted with the Worthey (1994) models. Worthey mod- 
els are currently available via the web at http : //astro . wsu . edu/worthey/dial/dial_a_model . html. 
F160W models were constructed by convolving the basic "vanilla" model spectral energy dis- 
tributions at each age and metallicity with the filter profile of the NIC2 F160W filter (G. 
Worthey, private communication). The default Worthey models use a Salpeter initial mass 
function and a helium fraction of Y = 0.228 + 2.7Z. There are differences in detail between 
the Worthey models and those presented in the previous sections. In particular, the Worthey 
models imply older ages for galaxies of a given color, and unreasonably large ages for the 
reddest galaxies. Note that the cahbration changes to the Freedman et al. (2001) Cepheid 
zero point without metallicity corrections makes the comparisons to the Worthey Mpieow 
models worse than they would have been using the previous Key Project calibration (Fer- 
rarese et al. 2002b, J2001). The comparison of the distance-independent (m/ — mpieow) 
models to the data yields ages and metallicities in excellent agreement with the conclusions 
of the previously-considered models. The redder galaxies have slightly sub-solar metallicities 
and old ages, while the bluer galaxies appear to have higher metallicities and significantly 
younger ages. It appears that the Worthey (1994) models may have an offset in absolute 
magnitude, but the relative brightnesses of the models at different wavelengths are reliable. 
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Fig. 3. — Absolute fluctuation magnitudes and fluctuation colors plotted with the Vazdekis 
stellar population models (Blakeslee et al. 2001), translated to the F160W fllter as described 
in the text. The symbols are the same as in Figure 2. 
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Fig. 4. — Absolute fluctuation magnitudes and fluctuation colors plotted with the Worthey 
(1994) models. The symbol deflnitions and uncertainties are the same as in Figure 2. Note 
that the triangular gap on the left is caused by the lack of published models for young, 
metal-poor populations, and is not a feature of the models. 
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4.4. Composite Stellar Population Models 

Fluctuation magnitudes are dominated by the youngest stars in a population. Statisti- 
cally, measured SBFs are the ratio of the second moment of the stellar luminosity function 
to the first, and therefore are very sensitive to the brightest stars. When comparison of SBF 
magnitudes to the models suggests the presence of young or intermediate- age populations 
(<5 Gyr), only a fraction of the stars would have to be young for the fluctuation amplitude 
to be high; the majority of the stars may be much older. 

Tonry, Ajhar, & Luppino (1990) were the flrst to attempt to use composite stellar 
population models to explain optical SBF observations, but the inadequacy of their model 
isochrones and colors led Worthey (1993) to question the vahdity of their results. More 
recently, Blakeslee et al. (2001) and Liu et al. (2002) have compared composite stellar 
population models to observational SBF data. Blakeslee et al. found that a 3-component 
model that includes a ~10% metal-poor component reproduces the observed /-band SBF 
slope and agrees with the IR SBF results presented here. 

In Figure 5 we compare our fluctuation magnitudes to the composite stellar population 
models of Liu et al. (2002) re-calculated for the F160W filter. In the three scenarios shown 
in Figure 5, 20% of the final mass formed six Gyr after the principal population. Gomposite 
populations with smaller mass fractions and smaller age differences would fall between the 
tracks shown in Figure 5; these models are meant to be representative and are shown for 
comparison with those published by Liu et al. (2002). The total time since formation is 
indicated in Figure 5 on each evolutionary track, thus the 7 Gyr point marks the magnitude 
and color of a population where 80% of the stars by mass are 7 Gyr old and the other 
20% formed 1 Gyr ago. The top line is for a burst with higher metaUicity ([Fe/H] = -1-0.4) 
than the original solar-met allicity population. The middle track is for a burst population 
of solar metaUicity, and the lowest track is for a burst with lower-than-solar metaUicity 
([Fe/H] = —0.7). In all three cases the majority of the stars (80% of the final mass) formed 
from solar-metallicity gas. 

Figure 5 shows that the bright fiuctuations in bluer galaxies may be the result of a recent 
burst of star formation within the last ~2 Gyr with a total mass fraction of ~20%. Liu et al. 
(2002) came to the same conclusion comparing A'-band SBF data to the composite models. 
The F160W data in Figure 5 suggest that the burst took place in previously-enriched gas. 
If recent episodes of star formation are triggered by mergers with dwarf galaxies with lower 
than solar metaUicity, the star formation must take place in gas previously enriched in the 
larger galaxy. Figure 5 suggests that as time goes on, bursts of star formation take place in 
increasingly metal-rich gas. 
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The bottom panel of Figure 5 shows the comparison of fluctuation color (m/ — mpieow) to 
the composite models. The (M/ — mFieow) models imply slightly lower absolute metalli cities 
for galaxies at the blue end of the color distribution than the Mpieow models do. None of 
the points in the lower panel fall much above the solar metallicity (m/ — mpieow) track, in 
contrast to the data in the top panel. The metallicities implied for the redder galaxies are 
the same in both panels. The metallicities implied by the (m/ — mpieow) models support 
the conclusion that more-recent star formation takes place in progressively higher-metallicity 
gas. 



4.5. Compciring Fluctuations to the H/3 Index 

The well-known agc-mctallicity degeneracy makes it difficult to use broad-band opti- 
cal colors to distinguish old, metal-poor populations from younger, metal- rich ones (e.g., 
Worthey 1994). One currently popular technique for determining ages and metallicities, 
particularly in young populations {t<5 Gyr), is to compare the age-sensitive H/3 absorp- 
tion line strength to a metal fine that is sensitive to metaUicity. IR SBFs also break the 
age-metallicity degeneracy and allow one to distinguish young, blue populations from older, 
metal-poor populations (Liu et al. 2000; Blakeslee et al. 2001). SBFs are dominated by 
the brightest stars in a galaxy, which in relatively young and intermediate-age populations 
(1 < t < 5 Gyr) arc evolved red giant and AGB stars. The populations we are considering are 
not actively forming stars, and are not so young as to contain massive luminous OB stars or 
emission nebulae. In contrast to SBFs, the H/3 absorption arises from main sequence stars 
near the main sequence turn-off. Because H/3 and SBF measurements arise from completely 
different groups of stars, and are measured using completely different techniques, the F160W 
SBF measurements presented here provide an important independent confirmation that H/3 
absorption is truly revealing differences in the ages of stellar populations. F160W fluctuation 
magnitudes and H/5 absorption are both sensitive probes of young stellar populations. The 
fact that the two techniques agree for individual galaxies in this sample is significant. 

In Figure 6 we compare the age-sensitive Mpieow to H/3 and use the Bruzual Sz Gharlot 
models of Liu et al. (2002) to infer relative ages and metallicities. We collected high-quality 
H/9 measurements on the Lick/IDS system from three recent studies (Kuntschner et al. 2000; 
Kuntschner 2000; Trager et al. 2000). These measurements include a sample of 14 galaxies 
in the Fornax cluster and 10 others, including some in the Local Group, Virgo, Leo clusters. 
The three studies measured H/3 values in apertures of different sizes. Kuntschner (2000) 
used a single aperture size for the Fornax cluster galaxies. Kuntschner et al. (2000) scaled 
their measurements by distance to a common physical scale. The Trager et al. (2000) data 




Fig. 5. — Three composite stellar population models (Liu et al. 2002) compared to the SBF 
measurements. In each case, the model population is composed of two bursts separated by 
6 Gyr, with the younger population containing 20% of the total mass. The total age of the 
composite population in Gyr is indicated by the numbers: 8 indicates the position on the 
line of a population where 80% of the stars by mass formed 8 Gyr ago and 20% formed 2 Gyr 
ago. The top (dashed) line indicates models with the second burst forming from super-solar 
metallicity gas ([Fe/H] = 0.4), the center (solid) line from solar metallicity gas ([Fe/H] = 0.0), 
and the lowest line (dotted) from sub-solar metallicity gas ([Fe/H] = —0.7). In all cases the 
older population is of solar metallicity. 
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are measured in an aperture scaled to the size of the galaxy (re/8). It is unnecessary to 
adjust the measurements to a common system since Kuntschner et al. (2000) found that 
H/3 does not change significantly with radius from the galaxy center. The Kuntschner and 
Trager studies apply a slightly different correction for H/? emission; we applied a small 
correction of 0.1 times [OIII]A5007 to the Trager et al. data to match the correction used 
in the Kuntschner papers. The H/5 measurements are listed in Table 1. Two galaxies were 
measured by both teams: H/3 measurements for NGC 3379 are presented in all three papers, 
and data for NGC 4472 are published in two. The H/3 values hsted in Table 1 are averages of 
all measurements, and the uncertainties were added in quadrature. The standard deviation 
of the different measurements was less than or equal to the total uncertainties in both cases. 

The trend seen in Figure 6 shows that the galaxies span a wide range in age, both as 
determined using H/3 and Mpieow- There is clearly a correlation between Mpigow ^-nd H/3. 
The stellar population models (Liu et al. 2002) have lines of constant age perpendicular to 
the correlation between Mpieow and H/9 seen in Figure 6, indicating that both are sensitive to 
young and intermediate- age populations (even though they sense very different components of 
the young populations). The fact that both techniques lead to the same conclusion that early- 
type galaxies often contain young stellar populations strengthens the conclusion considerably. 
The stellar population models become somewhat more compressed for ages greater than 8 
Gyr, when the A stars responsible for H/9 absorption have evolved off the main sequence, 
and the red giant and AGB stars responsible for the fluctuations are intrinsically fainter 
and more uniform in brightness (the model ages are best interpreted in a relative sense; 
observations that agree with the oldest models do not imply ages older than the Universe) . 

The absolute metallicities implied by Figure 6 are only marginally higher than those 
seen in the previous comparisons, and the trend with age is consistent. In the lower panel of 
Figure 6 the H/3 values are compared to the fluctuation color (m/ — mpieow), which is distance 
independent. The Mpieow and (m/ — mpieow) models disagree shghtly. The fluctuation color 
models suggest somewhat lower metallicities than the Mpieow models, in agreement with 
the other model comparisons. 

Both panels in Figure 6 show a relationship between IR fluctuations and H/3. The 
models suggest that both elliptical and SO galaxies span a very wide range in age and a 
relatively narrow range of metallicity. This result disagrees somewhat with the conclusion 
of Kuntschner (2000), who found that only the SO galaxies in Fornax have a wide spread in 
age. Kuntschner's limited sample of Fornax ellipticals appear uniformly old and span a range 
in metallicity. Our results are consistent with the conclusions of Trager et al. (2000), who 
found a wide age spread among a sample of ellipticals drawn from a variety of environments. 
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Fig. 6. — Absolute fluctuation magnitudes and fluctuation colors plotted vs. the H/3 absorp- 
tion line index (in A) for a subset of 24 galaxies (Kuntschner et al. 2000, and Kuntschner 
2000, and Trager et al. 2000). The majority of the galaxies are ellipticals in the Fornax 
cluster (see Table 1). H/?, like MpieoWj is sensitive to the age of the stellar population. 
Model predictions of both Mpieow and H/? are shown for the Liu et al. (2000, 2002) models, 
with ages and metalli cities indicated. The symbol definitions are the same as in Figure 2. 
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4.6. Common Ground 

In detail, there are significant differences between the different models in predicted ffuc- 
tuation magnitudes. Differences have also been noted at other wavelengths (Blalceslee et al. 
2001, Liu et al. 2002, Mei et al. 2001b). In general, however, the comparison of empiri- 
cal fluctuation magnitudes to the models yields a broadly consistent picture. Fluctuation 
amplitudes in galaxies redder than (y— 7)o>1.16 are all consistent with stellar populations 
older than approximately 8 Gyr and of slightly sub-solar metallicity. Galaxies on the bluer 
end of the color range all appear consistent with the youngest (2 to 5 Gyr) models in both 
absolute fluctuation magnitude and fluctuation color. Regardless of which model is consid- 
ered, these blue, lower-luminosity galaxies appear to be more metal-rich than the redder 
ellipticals. None of the bluer galaxies in this sample have fluctuation magnitudes that are 
consistent with old (>5 Gyr) single-burst stellar population models with low metallicities. 
The majority of stars in the blue galaxies may still be old, as the fluctuation magnitudes are 
dominated by young stars in a population. Composite stellar population models reinforce 
the conclusions of the single-age stellar population models. 

The conclusion that bluer early-type galaxies are younger and more metal-rich seems to 
contradict the observed relationship between mass and metallicity in giant ellipticals where 
less massive (bluer) systems retain less enriched gas, and should therefore have lower metal- 
licities than the massive eUipticals. The SBF magnitudes are dominated by the brightest 
and youngest stars in a galaxy, however, so our conclusion regarding metallicity may be 
consistent with the mass-metallicity relation if the bluer ellipticals have lower metallicities 
overall and some small fraction of their stars formed relatively recently from more metal-rich 
gas. The uncertainty in the assumed luminosity evolution of the AGB stars in the youngest 
population models also serves to minimize the importance of the apparent contradiction. A 
modest excess in the real AGB population over the models would result in enhanced SBF 
magnitudes that would appear the same as enhanced metallicities. 

The differences between models indicate that the details of the evolutionary tracks that 
are used in the models deserve further attention. Furthermore, the Padova evolutionary 
tracks are not computed beyond the onset of the thermally-pulsating stage of the AGB, and 
the modelers have employed a synthetic parameterization of the evolution of AGB stars to 
attempt to reproduce the luminosity function of realistic populations. The data presented 
here can provide feedback that should result in improved understanding of the luminosities 
and lifetimes of luminous AGB stars. It is also clear that real galaxies are not composed of 
coeval populations of stars of constant metallicity. 

The excellent agreement between the data and the models for both the Mpieow and 
the distance-independent {mj — mpieow) models suggests that the Cepheid distances are 
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probably reliable to better than 10%. For the current study we adopted the new Cepheid 
period-luminosity relations of Udalski et al. (1999), no metallicity corrections, and a distance 
modulus to the LMC of 18.50 mag. 

5. Galcixy Morphologies and SBF Magnitudes 

Figure 7 shows the absolute fluctuation magnitudes separated into bins by morphological 
type. Galaxy morphology classifications were taken from de Vaucouleurs et al. (1991, RC3). 
Absolute fluctuation magnitudes and the projected ellipticities of the galaxies on the sky are 
plotted in Figure 8 (RC3). The shapes of the symbols in Figure 8 are the same as the optical 
shapes of the galaxies measured well outside the NICMOS field of view. Morphological types 
and ellipticities from RC3 are listed in Table 1. 

As a group, lenticular SO galaxies are bluer than the elliptical population; nevertheless, 
the overlap between elliptical and SO galaxies covers almost the entire range in color covered 
by this sample. The bulges of spiral galaxies also span the full range of color and Mpieow 
spanned by the elliptical galaxies. Some dusty spiral bulges have fluctuation magnitudes that 
are significantly brighter than the early-type galaxies of the same color, most likely because 
the fluctuation amplitude is enhanced by the presence of clumpy dust. Some galaxies of 
all types included in our sample (T = —6 to -|-4) apparently contain bluer and brighter 
intermediate-age populations near their centers. While most of the galaxies are found in 
clusters or groups, the sample covers a range of environment from dense, compact clusters 
like Fornax, to lower-density groups and clusters (Leo, Virgo, and Ursa Major). Thirteen 
galaxies are not associated with any obvious group or cluster (as defined by Faber et al. 
1989). 

6. Summciry 

1. Using F160W SBF magnitudes to accurately determine extragalactic distances re- 
quires that the color of the galaxy be known. The best calibration fit is given in Equation 1. 
To achieve an accuracy of ~10% or better, the uncertainty in {V—I)o color should be less 
than ~0.035 mag. Applying the new calibration to the distant F160W SBF data of J2001 
would result in a value of the Hubble constant that is 10% larger than that published pre- 
viously; almost all of the difference (8%) is due to the improved Cepheid period-luminosity 
relation of Udalski et al. (1999). 

2. The reddest, most-massive ellipticals appear older (when compared to theoretical 
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Fig. 7. — Absolute fluctuation magnitudes plotted vs. color by galaxy type. T-types (RC3) 
of —4 to —6 are plotted as ellipticals, —3, —2, and —1 as lenticular (SO) galaxies, and types 
T>0 are shown as spirals. As in the previous figures, the largest symbols have the smallest 
uncertainties (median uncertainties are indicated in the upper right corner). 
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Fig. 8. — Absolute fluctuation magnitudes plotted vs. color, with symbol shapes that indi- 
cate the projected ellipticity on the sky (RC3). 
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stellar population models) than the less-luminous bluer galaxies. There are no galaxies in 
the blue half of our sample that have fainter fluctuation magnitudes consistent with very low 
metallicity ([Fe/H] < —0.7) and old age (t>5 Gyr) stellar population models. Comparison 
with stellar population models suggests that the youngest galaxies have somewhat higher 
metallicities than the older ellipticals. None of the galaxies have fluctuation magnitudes that 
are too bright to be accommodated by the young, metal-rich models. 

3. Composite stellar population models composed of ~20% younger stars by mass 
are also consistent with the brighter fluctuation magnitudes in bluer galaxies, provided the 
younger stars formed from gas of equal or higher metallicity than the gas that formed the 
older stellar population. 

4. The H/3 line index and fluctuation magnitude Mpieow are both sensitive to young and 
intermediate-age stellar populations with ages between one and a few Gyr. Age estimates 
from the two techniques are consistent. The agreement is significant because the two types 
of measurements make use of completely different techniques and are sensitive to different 
populations of stars. 

5. Comparison of fluctuation magnitudes and distance-independent fluctuation colors to 
the predictions of stellar population models provides an independent check on the distance 
scale calibration. If the models are correct, at least in the relative brightnesses between the 
/ and F160W bands, then the total systematic error in the Cepheid distance scale is possibly 
<0.1 mag. 

6. The SO galaxies and dust-free spiral bulges in this sample have fluctuation magnitudes 
that are indistinguishable from those of the eUiptical galaxies of the same color. Young and 
intermediate-age populations of relatively high metallicity must therefore be present in all 
galaxy types on the blue end of the color distribution. While there is a greater fraction of 
SO galaxies among the bluer galaxies, galaxies of all types are found across the entire color 
range. 

This study benefited greatly from NICMOS data collected as part of several programs, 
and we thank those who worked to acquire that data. In particular, we are grateful to those 
who helped ensure that the data would be appropriate for SBF analysis and assisted with the 
data reductions (D. Geislcr, J. Elias, J. R. Graham, and S. Chariot). Wc greatly appreciated 
the helpful comments of L. Ferrarese, T. Lauer, M. Postman, and R. Weymann. This study 
made use of data collected by the Optical SBF team (J. Tonry, J. Blakeslee, E. Ajhar, and A. 
Dressier) , and we thank them for providing color photometry and /-band SBF data. Finally, 
we wish to thank S. Chariot, A. Vazdekis, and G. Worthey for constructing the SBF stellar 
population models. 
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Table 1. F160W NIC2 Observational Data 
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26 
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33 
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68 
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07 
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32 


19 


8 
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36 
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21 
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04 
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1052 


207 


7886 


N4RW09 


02 


41 


04 


75 
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20 


7 
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-5 
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01 
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50 
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-4 
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-3 





40 


1 


50±0 


10 


NGC 


1339 


31 


7458 


N4B710 


03 


28 


06 


68 


-32 


17 


04 


3 


256 





057 


1392 


-4 
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34 


3 
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1 
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31 
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16 
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15 


56 


5 
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8 
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41 
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280 


558 


-3 
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NGC 


1427 


31 


7458 


N4B704 


03 


42 


19 
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23 


33 


8 
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-4 
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1 
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05 


NGC 
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09 
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06 


29 


2 


640 
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-5 
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048 


1364 


-5 
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10 
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640 
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37 
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3 


640 
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NGC 
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17 
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640 
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NGC 
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32 
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03 
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+2 
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6 


640 
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6 


640 
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NGC 
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28 


67 
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384 
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46 
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384 
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384 
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+3 
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37 
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640 





121 


1980 


-6 
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39 
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8 
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+2 





29 








NGC 


5193 




7453 


N4HD11 


13 


31 


53 


26 


-33 


14 


05 


2 


1920 
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640 
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09 
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06 
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640 
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NGC 7743 
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68 


09 


56 


03 
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640 





296 
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^Groups and clusters as defined by Faber et al. 1989 (Forna^=31, LeoI=57, Virgo=56); entries in parantheses are spirals not classified by Faber 
et al. but clearly belonging to a group or cluster. 

'^S-band extinction from Schlegel et al. 1998 

•^T-types and ellipticities from RC3; ellipticity e = (1-b/a) = (1-10~ ^""^ ^25) 
*^Kuntschner 2000, Kuntschner et al. 2000, and Trager et al. 2000 
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Table 2. Distances and F160W SBF Magnitudes 
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06 




new rL+Zj 


NGC 


404 


i 


Uo4±U 


Uii 


Z / 


41±0 


10 


21 


98±0 


09 


-5 


43±0 


13 
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Table 2 — Continued 
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31 


43±0 


29 


26 


61±0 


05 


-4 


82±0 


29 




7-SBF 


NGC 1553 


1 


159±0.016 


31 


18±0 


17 


26 


32±0 


07 


-4 


86±0 


18 




7-SBF 


NGC 3031 


1 


187±0.011 


27 


80±0 


26 


22 


99±0 


05 


-4 


81±0 


26 




7-SBF 








27 


75±0 


08 


22 


99±0 


05 


-4 


76±0 


09 




new PL 








27 


80±0 


08 


22 


99±0 


05 


-4 


81±0 


09 




new PL+Z 


NGC 3032 


1 


073±0.019 


31 


55±0 


28 


26 


05±0 


10 


-5 


50±0 


30 


D 


7-SBF 


NGC 3056 


1 


073±0.023 


30 


27±0 


25 


25 


48±0 


07 


-4 


79±0 


26 




7-SBF 


NGC 3351 


1 


225±0.014 


29 


85±0 


09 


25 


19±0 


07 


-4 


66±0 


11 


D 


new PL 








30 


00±0 


09 


25 


19±0 


07 


-4 


81±0 


11 




new PL+Z 


NGC 3368 


1 


145±0.015 


29 


92±0 


22 


24 


88±0 


09 


-5 


04±0 


24 


D 


7-SBF 








29 


97±0 


06 


24 


88±0 


09 


-5 


09±0 


11 




new PL 








30 


11±0 


06 


24 


88±0 


09 


-5 


23±0 


11 




new PL+Z 


NGC 3379 


1 


193±0.015 


29 


96±0 


11 


25 


26±0 


08 


-4 


70±0 


14 




7-SBF 


NGC 3384 


1 


151±0.018 


30 


16±0 


14 


25 


34±0 


17 


-4 


82±0 


22 




7-SBF 


NGC 3928 


1 


096±0.015 


31 


±0 


6 


25 


60±0 


09 


-5 


4 ±0 


7 


D 


7-SBF 


NGC 4143 


1 


181±0.015 


30 


85±0 


19 


26 


68±0 


06 


-4 


17±0 


20 




7-SBF 


NGC 4150 


1 


071±0.017 


30 


53±0 


24 


25 


60±0 


13 


-4 


93±0 


27 


D 


7-SBF 


NGC 4261 


1 


258±0.014 


32 


34±0 


19 


28 


04±0 


11 


-4 


30±0 


22 




7-SBF 


NGC 4291 


1 


175±0.017 


31 


93±0 


32 


27 


66±0 


08 


-4 


27±0 


33 




7-SBF 


NGC 4278 


1 


161±0.012 


30 


87±0 


20 


26 


38±0 


09 


-4 


49±0 


22 




7-SBF 


NGC 4406 


1 


167±0.008 


31 


01±0 


14 


26 


26±0 


06 


-4 


75±0 


15 




7-SBF 


NGC 4434 


1 


125±0.015 


31 


98±0 


17 


26 


78±0 


12 


-5 


20±0 


21 




7-SBF 


NGC 4458 


1 


140±0.011 


31 


02±0 


12 


26 


03±0 


05 


-4 


99±0 


13 




7-SBF 


NGC 4472 


1 


218±0.011 


30 


90±0 


10 


26 


26±0 


04 


-4 


64±0 


11 




7-SBF 


NGC 4527 


1 


23 ±0.03 


30 


53±0 


09 


25 


43±0 


07 


-5 


10±0 


11 


D 


new PL/G&S^ 








30 


61±0 


09 


25 


43±0 


07 


-5 


13±0 


11 




new PL/G&S+Z 


NGC 4536 


1 


20 ±0.07 


30 


80±0 


04 


25 


46±0 


12 


-5 


34±0 


13 


D 


new PL 








30 


87±0 


04 


25 


46±0 


12 


-5 


41±0 


13 




new PL+Z 
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Table 2 — Continued 



\JJ(XLcXA.y 




(V— 


(m-M) 
(mag) 


^Fieow 
(mag) 


-^Fieow 
(mag) 


i-/ Llot 




NGC 4565 


1 


128±0.027 


31 


05±0 


17 


25 


53±0 


08 


-5 


52±0 


19 


D 


7-SBF 


NGC 4589 


1 


180±0.015 


31 


55±0 


22 


27 


21±0 


08 


-4 


34±0 


23 




7-SBF 


NGC 4594 


1 


175±0.031 


29 


79±0 


18 


25 


06±0 


09 


-4 


73±0 


20 




/-SBF 


NGC 4636 


1 


233±0.012 


30 


67±0 


13 


26 


10±0 


08 


-4 


57±0 


15 




J-SBF 


NGC 4709 


1 


221±0.015 


32 


88±0 


17 


28 


51±0 


07 


-4 


37±0 


18 




HST J-SBF 


NGC 4725 


1 


209±0.023 


30 


45±0 


34 


25 


72±0 


10 


-4 


73±0 


35 




J-SBF 








30 


38±0 


06 


25 


72±0 


10 


-4 


66±0 


12 




new PL 








30 


46±0 


06 


25 


72±0 


10 


-4 


74±0 


12 




new PL+Z 


NGC 5193 


1 


208±0.015 


33 


35±0 


15 


28 


52±0 


06 


-4 


83±0 


16 


D 


HST J-SBF 


IC 4296 


1 


199±0.015 


33 


53±0 


16 


28 


82±0 


08 


-4 


71±0 


18 


D 


HST J-SBF 


NGC 5273 


1 


142±0.017 


30 


93±0 


26 


26 


13±0 


09 


-4 


80±0 


28 




J-SBF 


NGC 5845 


1 


124±0.012 


31 


91±0 


21 


27 


43±0 


07 


-4 


48±0 


22 




J-SBF 


NGC 7014 


1 


248±0.015 


33 


84±0 


15 


29 


06±0 


11 


-4 


78±0 


19 




HST J-SBF 


NGC 7280 


1 


105±0.009 


31 


77±0 


22 


26 


42±0 


09 


-5 


35±0 


24 




J-SBF 


NGC 7331 


1 


120±0.017 


30 


43±0 


17 


25 


43±0 


08 


-5 


00±0 


19 


D 


J-SBF 








30 


81±0 


09 


25 


43±0 


08 


-5 


38±0 


12 




new PL 








30 


84±0 


09 


25 


43±0 


08 


-5 


41±0 


12 




new PL+Z 


NGC 7457 


1 


104±0.009 


30 


45±0 


21 


25 


57±0 


09 


-4 


88±0 


23 




J-SBF 


NGC 7743 


1 


080±0.009 


31 


42±0 


17 


25 


99±0 


06 


-5 


43±0 


18 




J-SBF 



^J-SBF=Tonry et aL 2001, shifted by —0.16 mag to the cahbration derived using the 
Udalski et al. 1999 Cepheid period-luminosity relation; HST J-SBF=Jensen et al. 2001 
and Laucr ct al. 1998, also calibrated using the new period-luminosity relation; new 
PL=Frccdman et al. 2001; new PL/G&S=Gibson & Stetson 2001; +Z indicates distances 
derived using a metallicity correction of —0.2 mag dex~^. 

The Cepheid period-luminosity relation and metallicity corrections used by Gibson & 
Stetson were the same as those used by Frcedman et al. to determine the "new PL" and 
"new PL-I-Z" distances. The Gibson & Stetson distances assume a distance modulus for 
the LMC of 18.45 rather than the value of 18.50 adopted by Freedman et al. 2001. We 
have therefore added 0.05 mag to the Gibson & Stetson values. 
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